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Abstract 

We show that a generalized Landau theory for the smectic A and C phases exhibits a biax- 
iality induced AC tricritical point. Proximity to this tricritical point depends on the degree of 
orientational order in the system; for sufficiently large orientational order the AC transition is 3D 
XYdike, while for sufficiently small orientational order, it is either tricritical or 1st order. We 
investigate each of the three types of AC transitions near tricriticality and show that for each type 
of transition, small orientational order implies de Vries behavior in the layer spacing, an unusually 
small layer contraction. This result is consistent with, and can be understood in terms of, the 
"diffuse cone" model of de Vries. Additionally, we show that birefringence grows upon entry to 
the C phase. For a continuous transition, this growth is more rapid the closer the transition is to 
tricriticality. Our model also predicts the possibility of a nonmontonic temperature dependence of 
birefringence. 
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I. INTRODUCTION 

Since its discovery in the 1970's [1], the nature of the smectic A- smectic C transition has 
been a topic of great interest. Early work showed that many systems exhibit a continuous 
AC transition which could be described by a mean field model near tricriticality [2]. A 
tricritical point v with associated neighboring 2nd order and weakly 1st order transitions was 
later found jj], 4|. The origin of an AC tricritical point has been of significant interest, 
with two main mechanisms having been proposed. The first is the coupling of the tilt to 
biaxiality, which in chiral systems is related to the size of spontaneous polarization 0, 4]. 
The second is the width of the A phase 51]. Another mechanism, involving a coupling 
between tilt and smectic elasticity has also been proposed [fjj], but this seems less likely. 
Until now, a comprehensive theory that addresses the effect of biaxiality on the nature of 
the AC transition has not been produced. 

More recently, much attention has been given to de Vries materials, which exhibit an 
AC transition with an unusually small change in layer spacing and a significant increase in 
birefringence (associated with an increase in orientational order) upon entry to the C phase 
Some de Vries materials exhibit another unusual feature, namely a birefringence that 

n q 

varies nonmonotonically with temperature [8|, [9( ; in particular, the birefringence decreases 
as the AC transition is approached from within the A phase. De Vries materials generally 
seem to have unusually small orientational order and follow the phase sequence isotropic (/) 
- A - C. In several de Vries materials, the AC transition seems to occur close to tricriticality 



Separate theoretical models 
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ll| have been developed, each of which predicts the 



possibility of a continuous AC transition with the two main signatures of de Vries behavior: 
small layer contraction and increase in birefringence upon entry to the C phase. There are 
differences between the assumptions used in the models, the most significant of which is 
the treatment of the temperature dependence of the layering order parameter; the model of 
Gorkunov et al [11] does not take this into account while that of Saunders et al does 10]. 
Given the absence of a nematic phase in de Vries materials, incorporating the temperature 
variation of the layering order parameter is of crucial importance in the modeling of de Vries 
materials. It seems most likely that the I A transition in de Vries materials is primarily driven 
by the development of layering order, with orientational order being secondarily induced by 
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FIG. 1: Phase diagram in temperature (T) - concentration (c) space. For materials with excluded 
volume interactions, increasing the concentration would lead to an increase in the orientational 
order. The solid line represents the continuous AC boundary while the dashed line represents 
the 1st order AC boundary. These two boundaries meet at the tricrtical point: (T TC ,c TC ). The 
dotted line indicates the region in which the behavior in the C phase crosses over from Xy-like 
to tricritical. The region in which the behavior is Xl^-like shrinks to zero as the tricritical point 
is approached. Also shown as double ended arrows, are the three distinct classes of transitions (at 
fixed concentration): Xy-like, tricritical and 1st order. 

the layering order. This is consistent with the general observation |7| that de Vries materials 
have unusually strong layering order and unusually weak orientational order. Additionally, 



only by including temperature dependent layering, does one predict |lflj the unusual, yet 
experimentally observed [8J, |9j, possibility of a nonmonotonic temperature dependence of 
birefringence. 

Neither model considers the effect of biaxiality on the AC transition. The model of 
Gorkunov et al investigates the possibility of an AC transition that has signatures of tricrit- 
icality, but does not predict a tricritical point or the possibility of a 1st order AC transition. 



In this article, we present and analyze a new generalized nonchiral Landau theory, based 
on that developed in Ref. [loj |. which includes orientational, layering, tilt and biaxial order 
parameters. The model naturally produces a coupling between tilt and biaxiality and we 
show that this coupling leads to an AC tricritical point. We show that the effect of biaxiality 
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FIG. 2: A schematic showing the layer normal and optical axis. The layers are shown as dashed 
lines. The transition from the A to C phase occurs via a tilting, by angle 9, of the optical axis 
away from the layer normal. 

is stronger in systems with small orientational order, M , so that a tricritical point and 
associated neighboring 1st order transition can be accessed by systems with sufficiently 
small orientational order, Mq < Mtc- Here Mtc is the value of the orientational order at 
which the system exhibits a tricritical AC transition. This means that the two mechanisms 
that have been proposed as leading to tricriticality, the coupling of tilt to biaxiality and the 
width of the A phase, may in fact be two sides of the same coin. Systems with a narrow 
A phase, which are thus close to the I phase, will have small orientational order, which 
according to our model, leads to an enhanced effect of the biaxiality on the nature of the 
AC transition. For materials with excluded volume interactions, a decrease in orientational 
order could be achieved by decreasing concentration. 

Figure [1] shows the phase diagram for our model near the tricritical point in temperature 
(T) - concentration (c) space, along with the three different types of transitions: XV-like, 
tricritical and 1st order. In each case the transition from the A phase to the C phase implies a 
tilting of the optical axis away from the normal to the smectic layers by an angle 9, as shown 
schematically in Fig. [2j Our model gives the expected temperature dependence of 9 for each 
type of transition, as summarized in Fig. [31 For both the XV-like and tricritical transitions 
the growth of 9 with decreasing temperature is continuous, although with different scaling 
for each transition. It should be noted that here, and throughout the article, exponents 
are calculated within mean field theory, and do not include the effects of fluctuations. For 
example, it is known that when fluctuation effects are included in analysis of the 3D XY 
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transition, 9 scales like (1 — ^)^, with (3 ~ 0.35, whereas in mean field theory (3 = 0.5. The 
use of mean field theory is justified by the fact that virtually all continuous AC transitions 
are observed to be mean field like. 

For the 1st order transition the tilt angle 9 jumps discontinuously at the transition. 
Our model also leads to the expected [2j temperature dependence of specific heat c v near 
the continuous AC transition. This temperature dependence is shown in Fig. |H For an 
XY-like transition c v jumps by an amount Ac v as the system enters the C phase. If the 
transition becomes tricritical (M — > Mtc+, via decreasing concentration), the size of this 
jump diverges. Our model predicts that the divergence should scale like 

AC - * M -M TC ■ (1) 
For a 1st order AC transition there is an associated latent heat I. We show that if the 
transition becomes tricritical (M — > Mtc-, y i a increasing concentration) then the latent 
heat vanishes like 

I oc (M TC - M Q ) . (2) 



The model is also used to examine the behavior of the layer spacing and birefringence 
for the three possible transitions (XY-like, tricritical, 1st order). We show that, for all 
three types of transitions, an unusually small layer contraction can be directly attributed 
to unusually small orientational order, M . Specifically, we find that for any of the three 
possible types of transitions 

A d oc M (1 - cos(0)) « ^Motf 2 , (3) 

where the tilt angle 9 is small near a continuous or weakly 1st order transition. We define 
the layer contraction as A d = ((Iac ~ dc)/dAc, where (Iac an d dc are the values of the 
layer spacing in the A phase (right at the AC boundary) and in the C phase, respectively. 
Schematic plots of vs. 9 2 are shown in Fig. [5] for two types of systems: one " de Vries"- 
like and the other "conventional" . The "de Vries"-like system has small orientational order 
Mo <C 1 and thus has a small slope of vs. 9 2 , which corresponds to small layer contraction. 
The "conventional" system has strong orientational order M = 0(1), and thus has a larger 
slope, which corresponds to significant layer contraction. It should be noted that for a 1st 
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FIG. 3: The tilt angle 6 as a function of reduced temperature t=(l— jgj near the j4C transition 

temperature To, i.e., for t <C 1. Upon entry to the C phase the growth of the tilt angle scales like 

i i 
\t\ 2 for a mean field XY-like transition. For a tricritical transition it scales like \t\ 4 and is thus 

more rapid. For a 1st order transition there is a jump in the tilt angle upon entry to the C phase. 

order transition there will be a jump in the tilt angle 9 at the transition, and thus, the 
versus 9 2 line would not extend all the way to zero. 

This result of our rigorous theory complements the simple geometric diffuse cone ar- 
gument of de Vries 14j], which is shown in Fig. [6j The conventional, but oversimplified, 
relationship between layer contraction and tilt angle, A^ = (1 — cos(#)), is obtained geo- 
metrically by assuming a liquid crystal with perfect orientational order, as shown in Fig. 
EJa). However, it has long been known that the orientational order in liquid crystals is 
far from perfect. The schematic in Fig. [6](b) shows a more realistic arrangement of the 
molecules in the A phase. The molecular axes are tilted away from the optical axis, but 
in azimuthally random directions. One can see that the more the molecules are tilted, the 
smaller the orientational order in the A phase. The diffuse cone model argues that, upon 
entry to the C phase, the "pre-tilted" molecules do not need to tilt but rather need only to 
order azimuthally, thus leading to an unusually small layer contraction. Thus, the smaller 
the orientational order in the A phase, the more "pre-tilted" the molecules will be and the 
smaller the layer contraction will be. As shown in Eq. ([3]), our rigorous theoretical analysis 
predicts a small contraction for systems with small orientational order, which agrees with 
this geometric argument. It also correlates well with the general experimental observation 
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FIG. 4: The specific heat c v as a function of reduced temperature t = f 1 
AC transition temperature Tq, i.e., for t -C 1. As the transition is approached from C phase, the 



specific heat grows like c v oc [l — yf-J , where T m > Tq- This growth is cut off at T = Tc, 
where it reaches a maximum value, Ac v . If the transition becomes tricritical T m — > Tq and c v 
diverges at the transition. Note that the specific heat shown here only includes the contribution 
from the piece of the free energy density associated with the ordering as the system moves into the 
C phase. For a 1st order transition there will be a latent heat absorbed in going from the C phase 
to the A phase. 

j?! that de Vries materials have small orientational order. 

From Fig. MJd) one also expects a growth of orientational order, and hence birefringence 
An, as the system moves into the C phase. It is useful to define a fractional change in 
birefringence A An = A? ^^" AC , where Auac is the value of the birefringence in the A phase 
right at the AC boundary. Our model predicts that upon entry to the C phase, for any 
of the three types of transitions (XY-like, tricritical, 1st order), A An of a de Vries type 
material will grow according to A An oc 9 2 ■ While the dependence of A An on 9 is the same 
for all three types of transitions, its dependence on temperature is not the same because, as 
shown in Fig. [3l 9 scales differently with temperature for each type of transition. Thus, 
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FIG. 5: The layer contraction A^ = (cIac ~ dc)/dAc as a function of 9 2 near the AC transition. 
For any type of transition the contraction will scale like Mq9 2 . Thus, the slope of A^ versus 9 2 
is proportional to the orientational order Mq in the system. Near tricriticality, the orientational 
order is small and Mq <C 1 and so the contraction is also small. Also shown is the layer contraction 
for a system with strong orientational order Mq ~ 1, for which the contraction will be sizable. For 
a 1st order transition there will be a jump in the tilt angle 6 at the transition and thus, the A^ vs. 
9 2 line does not extend all the way to zero. 

The growth of Aa„ as a function of reduced temperature t = — 1 j is shown in Fig. 

[7J For an XY-like transition the growth will be linear oc \t\, while for a transition at 

i 

tricriticality it scales like oc \t\ 2 and is thus more rapid. For a 1st order transition there will 
be a jump in the tilt angle and thus an associated jump in Aa™, although near tricriticality, 
where the transition is only weakly 1st order, the jump will be small. 

Our model also predicts (for materials with excluded volume interactions) the possibility 
of birefringence that decreases as the AC transition is approached from the A phase, which as 
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discussed above, is an unusual feature that has been observed experimentally 
of the three types of transitions Aa w decreases linearly with temperature as the transition 
is approached from the A phase, as shown in Fig. UJ The decrease in birefringence is 
particularly unusual, as it indicates that the system is becoming less ordered (orientationally) 
as a lower symmetry (C) phase is approached. To the best of our knowledge, this is the first 
example of such a phenomenon. 

It should be emphasized that our analysis is only made tractable, and thus is only valid, 
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FIG. 6: (a) An oversimplified schematic showing the arrangement of molecules in the A phase, 
in which the orientational order is perfect. Such a model predicts that, as the system moves 
into the C phase, the layer spacing should contract according to A^ = (1 — cos(0)), where A^ = 
(d-AC — dc)/dAC- (b) A more realistic arrangement of the molecules in which the molecular axes 
are tilted away from the optical axis, but in azimuthally random directions. The more that the 
molecules are tilted, the smaller the orientational order. As the system moves into the C phase, 
the "pre-tilted" molecules do not need to tilt but rather need only to order azimuthally, thus 
leading to an unusually small layer contraction. Thus, the smaller the orientational order in the 
A phase, the more "pre-tilted" the molecules will be and the smaller the layer contraction will be, 
an interpretation consistent with our result, Eq. ([3]). The figure also shows that, as a result of the 
azimuthal ordering as the system moves into the C phase, it should become more orientationally 
ordered. 

in the limit of weak coupling between order parameters. This means that our results do 
not imply that all materials with small orientational order will have AC transitions close 
to tricriticality or will exhibit de Vries behavior. Similarly, not all materials exhibiting de 
Vries behavior must have AC transitions near tricriticality. In other words, the conclusions 
that our model leads us to are generic but not ubiquitous. The remainder of this article 
is organized as follows. In Section [Til we introduce our model and in Section II I II we locate 
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FIG. 7: The fractional change in birefringence A^„ = An ^^^ lc as a function of reduced temper- 
ature t = (l — Tjr^j near the AC transition temperature Tc, i.e., for t <C 1. For materials with 
excluded volume interactions, we expect the birefringence An, and thus Aa™, to decrease as the 
AC transition is approached from within the A phase. For all three types of transitions (XY-like, 
tricritical, 1st order) this decrease will scale linearly oc t with reduced temperature. Upon entry to 

the C phase the birefringence An, and thus A^ n , will grow. The growth is linear oc \t\ for a mean 

i 

field XY-like transition. For a tricritical transition the growth scales like oc \t\ 2 and is thus more 
rapid. For a 1st order transition there will be a jump in birefringence as the system enters the C 
phase. 



and analyze the biaxiality induced tricritical point. We then analyze the nature (XF-like, 
tricritical, 1st order) of the AC transition near this tricritical point in Section HVl In Section 
IVlwe examine the thermodynamic nature of each type of transition. Specifically, we calculate 
the specific and latent heats for the continuous and 1st order transitions, respectively. Lastly 
we study the behavior of the layer spacing and birefringence near the AC transition in Section 
ED We briefly summarize our results in Section IVIIi The Appendix includes details of the 
analysis from Section IVII 
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II. MODEL 



The starting point for our analysis is a generalized version of the free energy density 
introduced in Ref. jlQ| . which includes orient at ional, tilt (azimuthal), biaxial and layering 
order parameters. The complex layering order parameter ip is defined via the density p = po+ 
Re(?/>e Jcl ' r ) with po constant and q the layering wavevector, the arbitrary direction of which 
is taken to be z. The remaining order parameters are embodied in the usual second rank 
tensor orientational order parameter Q, which is most conveniently expressed as 

Qij = M[(-cos(a) + v / 3sin(a))ei i eij 
+(— cos(a) - v / 3sin(a))e 2 ie 2 j 

+2cos(a)e 3i e 3j ] , (5) 

where e 3 = c + y/1 — c 2 z is the average direction of the molecules' long axes, (i.e., the direc- 
tor). Here, in either smectic phase, z is normal to the plane of the layers. The projection, 
c, of the director onto the layers is the order parameter for the C phase. The other two 



principal axes of Q are given by ^ = z x c and e 2 = y/1 — c 2 c — cz. These unit eigenvec- 
tors are shown in Fig. [HJ The amount of orientational order is given by M oc x/Tr(Q 2 ), 
which is thus proportional to the birefringence. The degree of biaxiality is described by the 
parameter a. The A phase is untilted (c = 0) and uniaxial (a = 0), while the C phase 
is tilted (c ^ 0) and biaxial (a ^ 0). From Fig. [8] it can be seen that the angle 9, by 
which the optical axis tilts, can be related to c via c = sin(0). Taking both ip and Q to be 
spatially uniform allows the use of a Landau free energy density / = /q + + fg^, with 
the orientational (/q), layering (/^), and coupling (fq^) terms given by 

, t n Tr(Q 2 ) wTr(Q 3 ) , u n (Tr(Q 2 )) 2 



12 18 144 



(6) 



U = + \uM + \K{ q 2 - g 2 ) 2 |Vf, (7) 



fQ4> 



QiQj\1>\ 



-(a(q 2 )-b(q 2 M 2 )Q tJ +g(q 2 )Q lk Q jk 

h(q 2 ) s(q 2 ) 2 
H ^ — IkqiQkiQij 7 — (qkQlQkl) Qi 



(8) 



where the Einstein summation convention is implied and = qSi Z . As usual in Landau 
theory, the parameters t n and t s are monotonically increasing functions of temperature and 

11 



FIG. 8: The unit eigenvectors, ei, 62, S3 of the orientational order tensor Q. These are shown 
as solid arrows, with ei pointing into the page. Also shown, as a dotted arrow, is the layering 
direction z, which is normal to the plane of the layers. The eigenvector £3 corresponds to the 
average direction of the molecules' long axes. The order parameter, c, for the C phase is the 
projection of 63 onto the plane of the layers, and is shown as a dashed arrow. The angle 9, by 
which the optical axis tilts, is also shown. 

control the "bare" orientational and layering order parameters, M and ipo respectively. By 
"bare" we mean the values the order parameters would take on in the absence of the coupling 
term Jq^. Similarly, the constant go is the bare value of the layering wavevector. From Eq. 



(J7J) above, we immediately find \ipo\ = \/—t s /u s . The remaining parameters in fg and 
(w, u n , u s , K) are positive constants. 

The coupling piece of the free energy, fg^, includes the lowest order (in fields ip and 
Q) terms necessary to obtain an AC transition with tricriticality. The dependence on q 2 
of each of the coupling parameters, a, b, g, h and s, takes into account all other possible 
terms that have the same tensorial form, but with higher powers of q 2 , which is not an order 
parameter and is therefore not assumed to be small. For weak coupling, q « go we can 
Taylor expand each coupling parameter, e.g. a(g 2 ) ~ a + ai(g 2 — g^), where a = a(q 2 ), 



and a,i 



da 



For all but one of the couplings it is sufficient to use the zeroth order 

q z =Qo 



d(q 2 ) 

approximation, e.g. g(q 2 ) ~ go- It will be seen below that a±, the first order correction to 
do, is necessary for layer contraction at the AC transition. For notational convenience, we 
will, for the remainder of the article, write a(q 2 ) as a with the q 2 dependence implied. To 
render the analysis tractable, the coupling parameters are all assumed to be small and are 
treated perturbatively throughout. 

The relatively large number of parameters in / is inevitable given the fact that the theory 
incorporates four types of order, layer spacing and also allows for continuous, 1st order and 
tricritical AC transitions. Additionally, it will be shown that proximity to tricriticality and 
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the signatures of de Vries behavior can be interpreted simply in terms of the size of the 
orient at ional order. 



III. BIAXIALITY INDUCED AC TRICRITICAL POINT 

To investigate the nature of the AC transition, we expand the part of the free energy 
density involving orientational order, /q + Jq^ in powers of the biaxial and tilt order param- 
eters, a and c. This expansion is done near the continuous AC transition temperature Tq 
(i.e. for (T — Tc) /Tq <C 1) and to lowest order in M and ip- We find /q + /q* ~ fu + /coup- 
The piece fu only involves the orientational order parameter M and is given by 



From fu we immediately find the bare value of orientational order Mo(t n ) = (w + 
\Jw 2 — Au n t n ) /2u n . It is useful to write the orientational order as a combination of the 
bare value and a correction: M = Mq(1 + Am), where the correction A m is due to the 
coupling piece / COU p- The correction Am can be thought of as an augmentation of the bare 
orientational order M due to the presence of layering order. As discussed in Ref. [10], de 
Vries behavior is implied by a virtually athermal t n (and thus, an athermal M ), so that 
for a given material M can be thought of as a fixed quantity. This would correspond to 
almost perfect excluded volume short range repulsive molecular interactions. This means 
that the temperature variation in orientational order M is effectively due to its coupling to 
the temperature dependent layering, i.e. via Am- We assume and verify a posteriori that in 
the limit of weak coupling Am <C 1- Similarly, we express the wavevector as q 2 = q^l + A q ) 
and the layering order as \ip\ 2 = |^o| 2 (l + A^). The bare wavevector q is also taken to be 
athermal but the bare layering order parameter ipo is not. 

The coupling piece can be broken up into three pieces: / COU p = fuip + fc + fac- The piece 
/mi/; involves a coupling between layering and orientational order, that is non-zero in both 
A and C phases, and is given by 



I'm 



1 9 1 
-t n M 2 - -wM' 

2 3 




(9) 



fit* = q 2 \*P\ 2 M (-ar + g M - h q 2 M) , 



(10) 



where 



r = 1 



bo\ij\ 2 + (g + 2h q 2 )M 



(11) 



a 
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The piece f c involves the tilt (azimuthal) order parameter c and is given by 

fc = \r c c 2 + i Mc c 4 + X -v c c % . (12) 
The coefficients r c , u c , v c are given by 

r c = 3aq 2 \tfj\ 2 MT , (13) 

u c = 9h q 4 \^\ 2 M 2 , (14) 

v c = ^s q^\ 2 M 3 . (15) 

At the continuous AC transition the parameter r (and thus also r c ), changes sign. Close to 
the transition r oc (T — Tc)/Tc <C 1 and can be considered small. From Eq. (TTT]) we see 
that to lowest order in the corrections A.M,q,ip an d for athermal M , this transition, occurs 
due to layering order increasing as temperature decreases. The transition temperature Tq 



is defined via \ip (T c )\ = a/ (cl — (g + 2h qQ)M /b , or equivalently 

• (rp \ u s {a - {g + 2h q 2 )M ) 

% 

This continuous phase boundary is shown as a solid line in Fig. [9J the phase diagram in 
t s -M Q space. For a given material, decreasing the temperature would, in the phase diagram 
of Fig. [H correspond to moving horizontally from right to left. The size of the orientational 
order Mq should increase with concentration. Thus, the topology of the corresponding phase 
diagram, Fig. [fl in temperature-concentration space should essentially be the same as that 
shown in Fig. 

The coupling between tilt and biaxiality appears in the final piece 

fac = A a ac 2 + ^B a a 2 , (17) 



where, to lowest order in r, 



A a = 3 -^g oq 2 M 2 M 2 , (18) 
B Q = 3M 2 (wM - g q 2 \ij\ 2 ) . (19) 
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FIG. 9: The phase diagram in t s -Mo space near the tricritical point (t STC ,Mo TC ). The quantity Mq 
is a measure of how much bare orientational order the system possesses and for de Vries materials is 
effectively athermal. Increasing concentration should increase Mq. The quantity t s is a monotonic 
function of temperature so that for a given material, decreasing the temperature corresponds 
to moving horizontally from right to left. The topology of the corresponding phase diagram in 
temperature-concentration space should essentially be the same. The solid line represents the 
continuous AC boundary while the dashed line represents the 1st order AC boundary. These two 
boundaries meet at the tricritical point (t STC ,MQ TC ). The dotted line indicates the region in which 
the behavior crosses over from XYdike to tricritical. The region in which the behavior is XYdike 
shrinks to zero as the tricritical point is approached. The slopes of the 1st order and continuous 
AC boundaries are equal at the tricritical point. Also shown as double ended arrows, are the three 
distinct classes of transitions: XYdike, tricritical and 1st order. 

From Eq. (|17[) we see that biaxiality is induced by tilt order. Minimization gives 



Keeping in mind the weak coupling regime of our analysis, i.e. go <C 1, we see that the 
systems with small orientational order M will have large biaxial susceptibility. Thus, large 
biaxiality (and for chiral materials, an associated large spontaneous polarization) can be 
directly attributed to small orientational order. In fact, Eq. fT2Tj) predicts that the biaxial 



(20) 



where \ a can be thought of as a biaxial susceptibility and is given by 




(21) 
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susceptibility will be largest in systems that have a combination of weak orientational order 
(M) and strong layering order (|^|). It has been observed |7] that this combination may be 
common in de Vries materials. It should be noted that the expression for Xa is only valid for 
M > Ml = g q 2 \tfj\ 2 /w, below which terms we have neglected become important. However, 
we will see that the tri critical point we predict occurs at a value of M > Ml- 

The effect of the biaxiality on the AC transition is to renormalize the quartic coefficient 
in Eq. ffl2l) . giving 



For small biaxial susceptibility \a (corresponding to strong orientational order), the renor- 
malized quartic coefficient u' c > and the AC transition is continuous. For large Xa (corre- 
sponding to weak orientational order), u' c < and the transition is 1st order. The tricritical 
point occurs at r = u' c = 0, which, to lowest order in the corrections A^, corresponds to 
M = M TC with 



which is larger than Ml- For small coupling (ao, bo, go, h 1) the value of orientational 
order Mtc at tricriticality will also be small. In obtaining Eq. (1231) we have used Eq. 
( TTTj) at tricriticality to find I^OtcI 2 ~ a o/^o, an approximation that is valid for small M T q. 
Equivalently, t STC w —u s a /b . 

IV. AC TRANSITION NEAR THE TRICRITICAL POINT 

Having found the biaxiality induced tricritical point, we now investigate the nature of 
the AC transition in the vicinity of the tricritical point. We analyze both the continuous 
AC transition and the 1st order AC transition. 

A. Continous AC Transition Near Tricriticality 

For sufficiently large orientational order, M > Mtc, the renormalized quartic coefficient 
u' c > and the AC transition is continuous. As discussed in Section UTTl the phase boundary 
is defined via r = or equivalently t s = t s (Tc)- Upon entry to the C phase, r becomes 




(22) 




(23) 
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negative and, minimizing the effective f c (i.e. with u c — > u' c ) with respect to c we find that 
the tilt order parameter grows continuously with increasing \t\ like 



2h' Q 




9s q*M \ V (K 



2 



(24) 



where the effect of the coupling between biaxiality and tilt is incorporated via a renormalized 
h' , which by expanding \ a close to tricriticality (i.e. M ~ Mtc) can be shown to be 

Like u' c , h' changes sign at M = Mtc- It is straightforward to show that sufficiently close 
to the transition (|r| <C \ r *\), the dependence of c on r is effectively XF-like and that 
sufficiently far from the transition (|r| 3> |t*|) it is tricritical, i.e., 



O- ~~ ^/ 3h[^ 2 M (l T l) 2 



1 

~ ,, , , 1 

t\ > |r* 



(26) 



The crossover from XF-like to tricritical behavior occurs in the region r = 0(t*) where 
is the value of r where the c XY = c TC , 

M - ^ • (27) 

Near tricriticality where M is small, the corresponding t s * is given by i s * = t s (Tcr)(l + |r*|) 
and is shown as a dotted line in Fig. [91 The width of the region in which the behavior is 
XF-like shrinks to zero as the tricritical point is approached. Near the transition, the tilt 
angle 9 ~ c, and its scaling with temperature is shown in Fig. [3] for both an XlMike and a 
tricritical transition. Of course, the XY behavior of Eq. (126]) is that of a mean -field theory 
and incorporating fluctuation effects would yield c oc with f3 ~ 0.35. 



B. 1st Order AC Transition Near Tricriticality 



When the orientational order is small enough (M < Mtc) the quartic coefficient (u' c ) 
changes sign. The free energy now has two local minima, one at c = and another at 



2|^| 
9s q 2 M 



1 - 



4r 



(28) 
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The 1st order AC transition, and the jump from c = to c = c ut , occurs when the free 
energy at c lst becomes smaller than the free energy at c = 0. The location of the 1st 
order boundary can thus be obtained by finding where the two free energies are equal, or 
equivalently, where the difference Af between them is zero. To lowest order in corrections 
A.M, q ,ip this difference is just the effective f c (i.e. with u c — > u' c ) evaluated at c lat and is given 
by 




which when set to zero yields an expression for the location of the 1st order AC boundary 

3 



(30) 



Lst 16' 

This boundary is shown as a dashed line in Fig. [91 At the transition the tilt order parameter 



jumps from zero to a value c lstAc = ^\h' \/(3s q 2 M). Close to tricriticality, where the tran- 
sition is weakly lst order, c lgt is small and « 9. The corresponding temperature dependence 
of 9 is shown in Fig. [3l The size of the jump in c (and thus 9) goes to zero at the tricritical 
point, where h' Q — > 0_. 



V. THERMODYNAMIC NATURE OF THE AC TRANSITION NEAR TRICRIT- 
ICALITY 

We next investigate the thermodynamic nature of the AC transition near tricriticality. 
First we analyze the specific heat near the continuous transition and then the latent heat at 
the lst order transition. 



A. Specific heat near the continuous AC transition 

It is well established |2| that the specific heat will exhibit a jump at the continuous 
AC transition and that the thermodynamic signature of a continuous transition close to 
tricriticality is a divergence of this jump We obtain the specific heat for our model using 
c v = —T^ppf, where the prime indicates the use of the biaxiality renormalized u' c , as given by 
Eq. (|22|) . in f c . In using f' c instead of the full free energy density /, we are focussing on the 
contribution to the specific heat associated with the onset of ordering as the system moves 
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into the C phase. It is this contribution that is responsible for the specific heat jump. As 
discussed above, following Eq. ffl5|) . in a material with athermal Mo the transition from the 
A to C phase is driven by the layering order which increases with decreasing temperature. 
Near tricriticality, where the orientational order is small, the value of the layering order at 



the transition is l^oC^b)! ~ y a o/bo, and the dimensionless parameter r can be expressed as 

r^-mtLnJl-i), (31 ) 



\MTc)\ 2 ' \T C 

where we have Taylor expanded |^o(^)| near T = Tc and the dimensionless parameter 
7 C > is given by 7, 



Jb d\ip {T)\ 2 



expressed as 



. Using Eq. fl31|) . the specific heat can be 

T=T C 



T I %) 1? ' (32) 



In the A phase, where f' c = 0, the specific heat is zero. Using Eq. (|24p for c and Eq. fTTZj) 
(with u c — > u' c ) for f' c we can find the specific heat in the C phase. Thus we find 




T ) 2h' () 



4|t 1 1 * 1 V V l T *l 



r > 
r < 



(33) 



Close to tricriticality, where r* is small, the specific heat in the C phase near the transition 
is dominated by the first term. Substituting |r| = 7 C ^1 — ^7 J (valid in the C phase where 
T < Tc) into the first term, we find that c v scales like 



1 - jr) ■ (34) 

where T m = T c ^1 + > Tc. This scaling is shown in Fig. HI where it can be seen that 
specific heat grows as the AC transition is approached from the C phase. This growth is cut 
off at T = Tc (or equivalently r = 0), where it reaches a maximum value. This maximum 
value is the size of the specific heat jump at the AC transition and is found to be 

\T C J 2h' 

If the transition becomes tricritical then T m — > Tc and c v diverges at the transition. Equiv- 
alently, at tricriticality h' = and size of the jump Ac v diverges. Using Eq. (125!) we 
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Acy 



(M -M, c )- 



J j ^ M„ 

FIG. 10: The size of the specific heat jump Ac v as a function of the system's orientational order 
Mq. As Mo — > Mtc the transition becomes tricritical and the specific heat jump diverges. For 
systems with athermal Mo it should be experimentally possible to drive the system to tricriticality 
by varying the concentration. 

can relate a system's bare orientational order Mq to its proximity to tricriticality (where 
Mo = Mtc) which gives 

^«{wc- i T- (36) 

This relationship, shown in Fig. [TUJ allows us to see how the size of the jump in specific 
heat would diverge if the orientational order in the system could be tuned to approach Mtc- 
For systems with athermal M it should be experimentally possible to drive the system to 
tricriticality by varying the concentration. 



B. Latent heat at the 1st order AC transition 



For a 1st order AC transition there will be a latent heat absorbed in going from the C 
phase to the A phase. This latent heat vanishes when the transition becomes tricritical. We 
obtain the latent heat I for our model using I = —Tc^f: evaluated at the 1st order boundary, 
where for f c we use the expression given in Eq. ( J29l) . Using the relationship between r and 
T, as given in Eq. (IHT1) . we find 



I = lc-r- 
dr 



As the transition becomes tricritical h' — > 0_ and the latent heat vanishes. Relating the 
system's bare orientational order M to its proximity to tricriticality (where M = Mtc) 
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gives 



I oc 1 



M TC 



Mi 



o 



) 



(38) 



This relationship allows us to see how the latent heat would vanish if the orientational order 
in the system could be tuned to approach Mtc- For systems with athermal Mo it should be 
experimentally possible to drive the system to tricriticality, and the latent heat to zero, by 
varying the concentration. 

VI. BEHAVIOR OF THE LAYER SPACING AND BIREFRINGENCE NEAR THE 
AC TRANSITION 

We next analyze the behavior of the orientational order (which is proportional to the 
birefringence) and the layering wavevector (which is inversely proportional to layer spacing 
d) close to the AC transition. As discussed following Eq. (jUJ) above, for athermal M 
and go, the temperature variation of M — M (l + Am) and q 2 = q 2 (l + A q ) comes from 
the corrections Am and A q respectively. We thus seek the temperature dependence of the 
corrections Am i9 near the AC transition. Assuming, and verifying a posteriori, that the 
corrections are small, we Taylor expand the free energy to order (A^,,) 2 and minimize with 
respect to Am,q, keeping only terms to lowest order in coupling coefficients. This is done 
both within the A phase and within the C phase. Details of the analysis are given in the 
Appendix [A] 

A. Orientational order near the AC transition 

For the orientational order correction within the A phase we find 



where r is just the bare value of r, i.e., r evaluated at M = M , ip = ijj and q = q . To 
zeroth order in corrections A^ i9 , r = r . The quantity A° M = — 3(?oQ , oIV'o(7c')| 2 /7m < and 
for a continuous transition is just the value of the correction at the continuous AC boundary, 
i.e., where r = 0. At the 1st order AC boundary near tricriticality, at which r = T ut > 0, 
the correction is a little bit larger than A M [jj]]. Lastly, 7m = d 2 fui / dM 2 \ M=M . 




(39) 
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From Eq. (139|) we see that as the AC transition is approached from the A phase, i.e. 
as ro —>■ 0+, the correction Am a will decrease. For materials with sufficiently athermal Mo, 
this means that the orientational order will decrease as the transition is approached from 
above. Using the fact that birefringence An is proportional to orientational order M, the 
fractional change in birefringence A^ n = Ar ^^" AC (where the reader is reminded Auac is 
the value of the birefringence in the A phase right at the AC boundary) can be related to 
Am- It is straightforward to show that, to lowest order in Am, A^„ w Am — A° M . Thus, in 
the A phase A& n oc r will decrease as the transition is approached from above, as shown in 
7ig. (J7j). This is a feature that has been experimentally observed in some de Vries materials 
siull- We find this feature particularly interesting, as it is the first example that we know 
of in which the order of a phase decreases as a transition to a lower symmetry phase is 
approached. It should be noted that in materials with a sufficiently strongly temperature 
dependent t N , the growth of the "bare" (i.e., coupling-free) orientational order M (t n ) as T 
is lowered swamps the effects due to the correction term Am a - In this case, the orientational 
order would grow as the transition is approached from above. 

To find the correction near the transition within the C phase one must separately analyze 
the three distinct regions of the phase diagram, corresponding to XY, tricritical and 1st order 
behavior. As one might expect, the dependence of Am on r oc (T — Tc)/Tc -C 1 is different 
in each region. However, near tricriticality the dependence on the tilt order parameter c in 
each respective region (i.e. c XY , c TC and c lst ) is identical and is given by 

Kl |A 1 + If 1 + HM^A , (40) 



where A° M is equal to the value of the correction in the A phase right at the transition 15]. 
In each of the three regions the orientational order grows as one moves into the C phase, 
consistent with birefringence measurements of de Vries materials. Using the fact that the 
optical axis tilt angle 9 ~ c near the transition, we predict that the fractional change in 
birefringence will grow like A An oc 9 2 . It is important to note that while the dependence 
of the growth of A& n on 9 is the same in each of the three distinct regions of the phase 
diagram, the dependence on r is not. This is because the dependence of c (and thus 9) on 
ro differs in each of the three regions. For sufficiently large orientational order, away from 
the tricritical point c oc |t | i and the growth of A^ n near the continuous transition will 
scale like (Tq — T). For smaller orientational order, near the tricritical point c oc and 
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the growth of Aa„ will scale like (Tq — T) 2 . These scalings are shown in Fig. [7J Thus, our 
model predicts that for continuous transitions near tricriticality one will see a particularly 
rapid growth of birefringence as one moves into the C phase. For a 1st order transition there 
will be a jump in c and thus an associated jump in the birefringence. Close to the tricritical 
point, where the transition is weakly 1st order, this jump will be small. 



B. Layer spacing near the AC transition 

For the layering wavevector (which is inversely proportional to the layer spacing) within 
the A phase we find that 



A« = A ° + a j^r , (41) 



where A° = a\M§jK is value of the correction at the continuous AC boundary and the 



reader is reminded that a\ 



da 



. At the 1st AC boundary near tricriticality, at 



9 2 



dm 

which t = r l3t > 0, the correction is a little bit larger than A° [l6|. From the above 
equation we see that as the AC transition is approached, i.e. as tq — > 0+, the layering 
wavevector decreases. This corresponds to the layer spacing increasing, a feature which is 
generally observed experimentally. 

As with the orientational order, it is necessary to separately analyze the three distinct 
regions (XY, tricritical and 1st order) of the phase diagram to obtain the correction near 
the AC boundary in the C phase. Similarly, while the dependence of this correction on r 
differs within each region, the dependence on the respective tilt order parameter c in each 
region (i.e. c XY , c TC and c lst ) is identical. It is given by 

A, e =A; + »c>, (42) 

where A° is equal to the value of the correction in the A phase right at the transition 

n 

[1 61 ] and for a layer contraction (as opposed to dilation) to occur we have required a\ < 0. 
Using the above equation and the relationship between layer spacing (d) and wavevector 
(q = 2ft /d) we next seek the contraction in the layer spacing. This contraction is defined as 
Ad = (d AC — dc)/dACi where d A c and dc are the values of the layer spacing in the A phase 
(right at the AC boundary) and in the C phase respectively. We find that this contraction 
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is given by 

3|ai|M V. (43) 



2K 

Near the transition 9 ph c and the fractional contraction scales like 8 2 , as one would expect 
from the simple geometric argument discussed in the Introduction. However, our theory 
predicts that this fractional contraction is also proportional to the size of the orientational 
order, M ~ Mq. Thus, systems with unusually small orientational order will exhibit an 
unusually small layer contraction, as shown in Fig. [5j Given the fact that the tricritical 
point predicted by our model also occurs for small orientational order, it would not be 
surprising for some de Vries materials to exhibit AC transitions close to tricriticality. It 
should also be noted that for the 1st order transition, the contraction will be discontinuous, 
although the size of the discontinuity will nonetheless be proportional to the orientational 
order, which if small will make the contraction small. 



VII. SUMMARY 



In summary, we have shown that our generalized Landau theory exhibits a biaxiality 
induced AC tricritical point. The effect of the biaxiality is larger in systems with small 
orientational order, which would correspond to systems with narrow A phases. This means 
that the two mechanisms that have been proposed as leading to tricriticality in a system, 
the coupling of tilt to biaxiality and the width of the A phase, can both be attributed to the 
system possessing sufficiently small orientational order. For materials with excluded volume 
interactions, one could reduce the orientational order, and thus tricritical point, by 

reducing concentration. We have shown that the optical tilt, specific heat and latent heat 
all exhibit the expected behavior near tricriticality. In addition, we have explored the effect 
of proximity to tricriticality on these quantities, and we have quantified the effect in terms 
the degree of orientational order in the system. 

We have also analyzed the behavior of the birefringence (via the orientational order) 
and the layer spacing (via the wavevector) for each of the three possible types of transitions 
(XY-like, tricritical and 1st order) near tricriticality. For de Vries material the birefringence 
has been shown to increase upon entry to the C phase and for a continuous transition this 
increase is more rapid the closer the transition is to tricriticality. It was also shown that for 
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materials with excluded volume interactions, birefringence will decrease as the AC transition 
is approached from the A phase, implying a non-monotonic temperature dependence of 
birefringence, a very unusual feature. We have used our model to obtain a relationship 
between the layer contraction and the tilt of the optical axis as a system moves into the 
C phase, for any of the three possible types of transitions. This relationship predicts that 
systems with small orientational order in the A phase will exhibit a corresponding small 
layer contraction. Our result correlates well with the diffuse cone geometric argument of de 
Vries. 

Our future work in this area will involve further generalizing our model to include chirality. 
Having done so, we will analyze the electroclinic effect in materials near the AC transition. 
Of particular interest will be how the size of electro-optical response depends on orientational 
order and proximity to a tricritical point. 
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APPENDIX A: CORRECTIONS TO THE BARE ORIENTATIONAL ORDER 
AND TO THE BARE LAYERING WAVEVECTOR 

In this Appendix we outline the procedure by which we obtain the corrections, Am and 
A q , to the bare orientational order and to the bare layering wavevector, respectively. This is 
done near the AC boundary for both the A phase and the C phase. Near the AC boundary 
within the C phase, we analyze separately the three regions of interest (XF-like, tricritical 
and 1st order). 

1. Correction to the bare orientational order 

In this section we find the correction A M to the bare orientational order M , where A M 
is defined via the full orientational order M = M (l + Am)- This is done by expanding the 
free energy to order (A M ) 2 in the phase of interest and then finding the A M that minimizes 
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the free energy. 



a. Correction in the A phase 

We begin our analysis of the correction in the A phase by expanding /m, given by Eq. 



fu ~ f M {M ) + l - lM Ml^ M f , (Al) 



where 7m = d?fM/dM 
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\M=M - 



In both the A and C phases, a non-zero Am is due to the coupling parts of the free 
energy In the A phase only the piece /m^, given by Eq. ( flOl) . is non-zero. Expanding /mi/>, 
which requires the expansion of r, yields 

fit* « /W + 9oW 2 Mo (3c/ M - a r ) A M , (A2) 

where /m^ and r are the bare values of fu* and r, i.e. evaluated at M = M , ip = ip and 
q = go- We have ignored order (Am) 2 terms, which are higher order in the coupling than 
the (Am) 2 term in Eq. ( 1 All) and are thus subdominant. Minimizing /m + /mv> with respect 
to Am gives 

A Ma = — 77 {-3g M + a T ) , (A3) 

Mo7m 

where we have replaced ~ V'o(Tc') near the AC transition. The above expression can be 
rearranged to give Eq. fl3"§|) . From the above expression we see that the correction Am is 
on the order of the coupling parameters, and go, and is thus small as was assumed in 
expanding the free energy. 



b. Correction in the C phase 

In finding the corrections in the C phase near the AC boundary we first follow the same 
procedure as for the A phase, namely the expansion of fu and fu* as given by Eqs. (1A1[) 
and 1IA2P above. We must also expand the piece of coupling, f' c , that is non-zero in the C 
phase. The prime indicates the use of the biaxiality renormalized u' c , as given by Eq. 11221) . 
in f c , which is given by Eq. ([121) . For each separate region of interest (XY, tricritical and 
1st order) we use the appropriate expression for c in f' c . 
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In the XF-like region we find 



ri \ib\ 2 a 2 T 2 



fcxY 4< 4h' " (A4) 

Expanding r and h' in powers of Am, keeping terms to lowest order in tq and coupling 
coefficients gives 

4. « /cxy n + + 2h %) Am , (A5) 



where f' CXYo and /iq are the bare values of f' CXY and h' Q . 
Minimizing f M + f M f + f' CXY with respect to A M gives 

where, in neglecting the To dependent contribution from fu^-, we have used the fact that 
close to tricriticality h' 00 /ho <C 1. Using the bare version of c = c XY as given by Eq. fl26|) 
this expression can be rearranged to give Eq. PU|) . 

For the tricritical region where, u' c is effectively zero, one must use f' c evaluated at c = c TC 
which yields 



f, 1 ~r 3 c 2|^| 2 -a 3 r3 



Expanding r in powers of A M while keeping terms to lowest order in r and coupling coef- 
ficients gives 



4c « 4o - \MTc)\ 2 M ^j ( go + 2h q 2 ) A M , (A8) 

where /c TCo * s ^ ne ^ are vame °f f' CTC - 

Minimizing f M + /mv> + f' CTC with respect to A M gives 

where, in neglecting the r dependent contribution from /m^, we have used the fact that 
S> To close to tricriticality, i.e. where r* <C 1. Using the bare version of c = c TC as given 
by Eq. (1261) this expression can be rearranged to give Eq. (HOI) . 

Lastly we obtain the correction in Am in the C phase (where h' < 0) near the 1st order 
AC boundary. We do this by expanding f' c near the first order AC boundary, the expression 
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for which is given by Eq. (|29|) . Expanding r, h' and r* (which depends on h' Q ) in powers of 
Am while keeping terms to lowest order in tq and coupling coefficients gives 



fL - /;«„ - ^« (l + ^1 - ( 9 o + 2M3 A M , (MO) 

where f ' is the bare value of f ' . 

Minimizing / M + / M ^, + with respect to A M gives 



- (- 3 ^° + H i 1 + V 1 " £i) {s " + H ' (A11) 

where r* is the bare value of r* and, in neglecting the r dependent contribution from Jm^, 
we have used the fact that close to tricriticality h' 00 /ho <C 1. Using the bare version of 
c = c lgt as given by Eq. (128]) this expression can be rearranged to give Eq. (HO]) . 



2. Correction to the bare wavevector 



In this section we find the correction A q to the bare wavevector q , where A q is defined 
via the full wavevector q 2 = q 2 (l + A q ). As with the orientational order, this is done by 
expanding the free energy to order (A q ) 2 in the phase of interest and then finding the A q 
that minimizes the free energy. 



a. Correction in the A phase 

We begin our expansion of the free energy in powers of A q by expanding fy, given by 
Eq. ©, 

U * ±K\ip \Y A 2 q . (A12) 

In both the A and C phases, a non-zero A g is due to the coupling parts of the free energy. 
In the A phase only the piece /mi/>, given by Eq. (fit)}) , is non-zero. Expanding fu^ yields 

/W ~ /W - ql\^o\ 2 M (a x ql + a r ) A q , (A13) 

where we have used the fact that M is small near tricriticality. We have ignored order (A q ) 2 
terms, which are higher order in the coupling than the {A q ) 2 term in Eq. (1A12I) and are 
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thus subdominant. Minimizing Jm + fMip with respect to A q gives 

A ia = -77^2 ( ai <?o + a o T o) ■ ( A14 ) 

The above expression can be rearranged to give Eq. (j4"Tl) . From the above expression we 
see that the correction A 9 is on the order of the coupling parameters, a and ai, and is thus 
small as was assumed in expanding the free energy. 



b. Correction in the C phase 

In finding the corrections in the C phase near the AC boundary we follow the same 
procedure as for the orientational order. To obtain the correction within the XF-like region 
we use f' Cxy as given by Eq. (1A4I) . Expanding r and h' in powers of A q , keeping terms to 
lowest order in To and coupling coefficients gives 

f > „ f \MTc)\ 2 aiq 2 aoTo 

Jc XY ~ Jc XYo 7fT, > l AiD J 

where we have used the fact that M is small near tricriticality. 
Minimizing f M + f Mi> + f' CxY with respect to A q gives 

v^H-f). (AM) 

where, in neglecting the To dependent contribution from fu^-, we have used the fact that 
close to tricriticality h' 00 /ho 1. Using the bare version of c = c XY as given by Eq. (f26|) 
this expression can be rearranged to give Eq. (|42p . 

For the tricritical region we use f' CTC as given by Eq. flA70 . Expanding a and r in powers 
of A q while keeping terms to lowest order in tq and coupling coefficients gives 



fcrc « fU + \MTc)\ 2 q 2 a n l ^A q , (A17) 



3s 

where we have used the fact that M is small near tricriticality. 
Minimizing f M + f M ^ + f' CTC with respect to A q gives 
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where, in neglecting the To dependent contribution from fu^-, we have used the fact that 
^Jto 3> t close to tricriticality, i.e. where r* <C 1. Using the bare version of c = c TC as given 
by Eq. (125]) this expression can be rearranged to give Eq. f!42p . 

We conclude by obtaining the correction in A q in the C phase (where h' Q < 0) near the 
1st order AC boundary. We do this by expanding f' c near the first order AC boundary, the 
expression for which is given by Eq. (1291) . Expanding r and h' in powers of A q , keeping 
terms to lowest order in tq and coupling coefficients gives 

/;..-/;..„ (i+^g)A„ (A 19) 

where we have used the fact that M is small near tricriticality. 
Minimizing f M + f Ml p + /' with respect to A q gives 



Aifc - = ik [ Moql ~ 1 1 + (/l " ^ 1 1 • (A20) 




where, in neglecting the To dependent contribution from fu^-, we have used the fact that 
close to tricriticality h' 00 /ho 1. Using the bare version of c = c lst as given by Eq. 
this expression can be rearranged to give Eq. 
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